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Background: Nuclear fission is a complex large-amplitude collective decay mode in heavy nuclei. Microscopic 
density functional studies of fission have previously concentrated on adiabatic approaches based on constrained 
static calculations ignoring dynamical excitations of the fissioning nucleus and the daughter products. 

Purpose: We explore the ability of dynamic mean-field methods to describe fast fission processes beyond the 
fission barrier, using the nuclide 240 Pu as an example. 

Methods: Time-dependent Hartree-Fock calculations based on the Skyrme interaction are used to calculate 
nonadiabatic fission paths, beginning from static constrained Hartree-Fock calculations. The properties of the 
dynamic states are interpreted in terms of the nature of their collective motion. Fission product properties are 
compared to data. 

Results: Parent nuclei constrained to begin dynamic evolution with a deformation less than the fission barrier 
exhibit giant-resonance-type behaviour. Those beginning just beyond the barrier explore large amplitude motion 
but do not fission, whereas those beginning beyond the two-fragment pathway crossing fission to final states which 
differ according to the exact initial deformation. 

Conclusions: Time-dependent Hartree-Fock is able to give a good qualitative and quantitative description of 
fast fission, provided one begins from a sufficiently deformed state. 


I. INTRODUCTION 

Studies of nuclear fission have been ongoing since the 
discovery of the process in 1938 by Hahn and Strassmann 
[1]. The actinide nuclide 240 Pu has long been a case of 
interest, as spontaneous fission presents itself as a decay 
mechanism with significant probability, relative to other 
isotopes in the actinide region. This allows for quan¬ 
titative comparisons between spontaneous and induced 
fission [2-5] . Experimentally, fission can be induced by a 
variety of techniques, including neutron-induced fission, 
fission induced by more complex projectiles, and photo¬ 
fission [2, 6]. Recent experimental campaigns have inves¬ 
tigated /3-delayed fission [7]. 

Theoretically, microscopic studies have focused upon 
the role of the quadrupole degree of freedom in forming 
the fission pathway, as exemplified by constrained mean- 
field calculations [8, 9]. The typical observed behaviour 
in actinide nuclei for the binding energy as a function of 
increasing quadrupole deformation is to follow a multi¬ 
humped pathway (see Fig. 1). When considering the po¬ 
tential energy surface (PES), starting from the ground 
state, then increasing the quadrupole deformation will 
result in a first fission barrier. By increasing the de¬ 
formation further, a secondary minimum, corresponding 
to an isomer, is found. Beyond this minimum, a sec¬ 
ond fission barrier is encountered, and past this barrier, 
the general consensus is that it becomes more energet¬ 
ically favourable for the nucleus to fission. The ener¬ 
gies Ea, E b and Ejj presented in Fig. 1 correspond to 
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FIG. 1. (Color online) Schematic of typical potential energy 
surface obtained when increasing the elongation (which cor¬ 
responds to quadrupole deformation) of an actinide nucleus. 
The energies Ea , Eg and Eji relate to properties of the fis¬ 
sion barriers, and correspond to the definitions in Ref. [10]. 


those defined in Ref. [10]: the energy difference between 
the ground state and the peak of the first fission barrier, 
the difference between the ground state and the peak of 
the second fission barrier, and the difference between the 
ground state and fission isomer, respectively. In some ex¬ 
otic cases, triple-humped potential surfaces are expected 
[11-14]. Although the multi-humped behaviour of the 
energy surface cannot be measured directly, experimen¬ 
tal evidence points towards this characteristic structure 
[13, 15], 

For fission studies, the quadrupole degree of freedom 
is of capital importance, as it describes the elongation 
of the nucleus [16]. Additionally, as many nuclei are ob- 
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served to fission with asymmetric mass distributions, the 
octupole degree of freedom is vital to describe any mass 
reflection asymmetry. Modern density functional theory 
(DFT) solvers are able to perform symmetry-unrestricted 
calculations which allow, in principle, any and multiple 
degrees of freedom to be explored [17]. As well as con¬ 
sidering deformation degrees of freedom, constraints can 
be imposed from alternate perspectives to study fission 
in static calculations. Some studies, for instance, assume 
a symmetric fission fragment path [18] or a constrained 
multi-configurational static solution before investigating 
time evolution [19]. In fact, time-dependent generator 
coordinate method calculations provide a generalization 
of the approach presented here to take into account the 
evolution of collective coordinates [20] . 

The approach of calculating the PES to describe fis¬ 
sion, regardless of the number of particular degrees of 
freedom constrained, is limited to producing a series of 
static solutions which attempt to describe a dynamic 
process, resulting in an effectively adiabatic approxi¬ 
mation. Shape-constrained DFT calculations produce 
Slater determinants which contain no internal excita¬ 
tions. Some attempts have been made to account for 
finite-temperature effects as sources of dissipation [21]. 
How such effects would couple with the nonadiabatic 
time-dependent Hartree-Fock (TDHF) approach is an 
open question that deserves further study [22]. Time- 
dependent techniques coupled with constrained calcula¬ 
tions may therefore yield new, insightful results, because 
they describe the dynamics of a fissioning system. Time- 
dependent Hartree-Fock [23] presents itself as a candidate 
method, as it is able to temporally evolve Slater determi¬ 
nants which begin as a solution to the constrained static 
Hartree-Fock (HF) equations. 

The calculational basis for this study is the so-called 
TDHF technique. We note, however, that nowadays this 
should probably be called “time-dependent density func¬ 
tional” , because it is based upon a nuclear energy density 
functional. The slightly less correct TDHF moniker has 
stuck within the nuclear physics literature, though, and 
we keep this convention in the present work. 

TDHF is the basic lowest-order microscopic dynami¬ 
cal mean-field theory, first proposed by Dirac [23] and 
later on applied to more or less realistic nuclear systems 
in the 1970s—1980s [24-28]. A practical implementation 
involves beginning from an energy density functional and 
using the variational principle to obtain Hartree-Fock- 
like equations for the static initial state. The time evo¬ 
lution equations, based on the same energy density func¬ 
tional, are then run on this initial state. We use the 
Skyrme energy density functional, depending on the lo¬ 
cal densities and currents, 

B = E s ky{p, «7, t, s, j , £) (1) 

where p is the particle density, J is the vector part of the 
spin-current tensor, r is the kinetic density, s is the spin 
density, j is the particle current and £ is the pairing den¬ 
sity [29]. These densities and currents include time-odd 


fields (a and j), which are only active in the dynamic 
part of the calculation. Only those time-odd fields that 
couple to the necessary time-even fields through Galilean 
invariance are included in the present calculation. Full 
details of the practical aspects of solving the static and 
dynamic HF equations, along with the explicit details of 
the density functional are contained in the documenta¬ 
tion of the sky3d code [30]. 

There have been several historic attempts to describe 
fission dynamics using TDHF. The pioneering attempt 
by Negele et al. [26] was followed by a series of studies 
that suffered from computational power limitations and 
were hindered by axial symmetry and restricted forms 
of the nucleon-nucleon interaction [31, 32]. With a new 
generation of TDHF solvers able to perform symmetry- 
unrestricted, three-dimensional calculations [30, 33, 34], 
modern TDHF studies have begun to take a renewed in¬ 
terest in fission [18, 35, 36]. Here we use a modified 
version of the recently published code SKY3d [30], with 
shape constraints explicitly included. We then solve the 
HF + BCS equations in a three-dimensional Cartesian 
basis and evolve the calculated states using TDHF. 

Spontaneous fission cannot be accessed directly within 
TDHF. To reach a fissioned configuration from the 
ground or isomeric state, the nucleus must tunnel 
through the barriers in the PES. While TDHF allows a 
quantum mechanical description of single-particle wave 
functions, the collective motion is semiclassical, hence 
forbidding tunnelling in collective coordinates [37]. In 
contrast, TDHF is suitable for exploring the dynamics 
of induced fission. The potential challenging issue is 
the incorporation of the fissioning mechanism within the 
TDHF framework. Here we follow the strategy of find¬ 
ing constrained Hartree-Fock (CHF) states and use them 
as initial conditions in a TDHF calculation. Our aim 
is to investigate how the underlying deformed structure 
pushes the parent nucleus towards fissioning paths, if it 
does so at all. 

We also run TDHF simulations from initial states that 
lie past the fission barrier. At that stage, one is some¬ 
how mimicking spontaneous decay in the sense that these 
represent the states right after tunnelling. The complex 
quantum dynamics that occur in the lead-up to the point 
at which we let the TDHF calculations take over will 
clearly populate a range of different configurations, lead¬ 
ing to a spread of fission products. Our dynamical cal¬ 
culations only take a handful of initial states on a single 
constrained quadrupole PES, as a proxy for the quan¬ 
tum nature of the “tunnelling” process. A more realistic 
way of sampling may be to find multiple PES solutions 
at the same energy and start from those. We note that 
our different initial configurations correspond to slightly 
different total energy content in the system. 

The paper is organised as follows. In Sec. II, the static 
constrained starting points are discussed. Section III 
explores the dynamics of the fissioning nucleus, while 
Sec. IV looks in more detail at the dynamics of the fission 
fragments. Some concluding remarks are given in Sec. V. 
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II. STATIC CONFIGURATIONS IN 240 PU 

We begin the investigation of fission by examining the 
PES for a nucleus of interest. The reproduction of the 
double-humped fission barriers of actinide nuclei are of¬ 
ten used as a benchmark test for nuclear models [38, 39]. 
Owing to the wealth of data available from experimental 
[2-5] and theoretical [9, 10, 40, 41] studies, 240 Pu presents 
itself as a strong candidate for a benchmark test of TDHF 
to investigate induced nuclear fission [39]. We note, in 
particular, that configuration mixing plays a relatively 
small effect in this isotope [42]. The static quadrupole- 
constrained PES is firstly calculated to provide a selec¬ 
tion of initial states for time evolution. These are ob¬ 
tained by performing an energy minimisation with re¬ 
spect to constraints imposed upon the quadrupole shape 
degree of freedom. 

Modern symmetry-unrestricted DFT solvers have ex¬ 
tended the PES for multiple constraints, for example si¬ 
multaneously constraining the quadrupole and octupole 
degrees of freedom to explore two-dimensional deforma¬ 
tion surfaces [43]. The approach of calculating a multi¬ 
dimensional PES to describe fission within a microscopic 
framework has enjoyed much recent attention [14, 21, 43- 
45]. Alternatively, one can explore fission pathways via 
shell-corrected macroscopic liquid drop models. This 
technique has been applied to perform exhaustive topo¬ 
graphical surveys of deformation space to deduce fission 
properties of static configurations [10, 46-48]. However, 
this approach, regardless of the number of dimensions, 
is limited to producing a series of static solutions to de¬ 
scribe a dynamic process, thus being equivalent to adia¬ 
batic motion. Time-dependent simulations describe the 
fissioning system, while allowing for internal excitations. 
This is important for the final, fast stages of fission, which 
we explore in the present work. 

To build shape-constrained ground states, we begin 
from an arbitrarily deformed state and use the aug¬ 
mented Lagrangian method [43, 49] to constrain the 
quadrupole deformation. Our purpose is not to pursue 
an in-depth investigation of multiple shape-constraints 
in HF calculations, but rather to use the technique to 
produce initial states to then investigate their time evo¬ 
lution. All other degrees of freedom are assumed to settle 
into the configuration of minimum energy [37]. 

For the results presented here, where only one con¬ 
straint is applied upon the quadrupole degree of free¬ 
dom, a masking procedure has been adopted. This limits 
the space which the nuclear wave functions can explore, 
allowing a single fission pathway to be explored where 
the nuclear shape gradually evolves, rather than abruptly 
jumping between competing energy minima. Details of 
the masking procedure are discussed in Ref. [50]. Our 
emphasis is on the fast fission dynamics beyond the bar¬ 
rier, but it may yet be fruitful to use other constraints to 
generate more starting points. 

The ground-state and CHF calculations are performed 
using the SkM* effective interaction. The fission bar¬ 


TABLE I. Summary of ground-state and isomer properties 
of 240 Pu, calculated using SKy3d with the SkM* Skyrme ef¬ 
fective interaction. Further details of the calculations are in¬ 
cluded in the text. Both the ground state and isomer are 
prolate deformed and axially symmetric. 


Nucleus 

Binding Energy 
[MeV] 

rms Radius 

M 

P 20 

Aso 

A*0 

240 Pu 

240 Pu* 

-1781.95 

-1778.91 

5.941 

6.418 

0.280 

0.682 

0.000 

0.000 

0.255 

0.547 


rier properties of 240 Pu were considered when fitting 
the SkM* effective interaction [51]. We perform our 
static calculations in a regularly spaced Cartesian grid 
of 40 x 40 x 40 points, ranging from —19.5 to 19.5 frn 
in the x, y and z directions. This rather coarse grid 
gives surprisingly good results, with, e.g., binding ener¬ 
gies typically differing by parts in around 10 4 relative to 
much finer grids [30] . BCS pairing is included within the 
static calculation, using the Volume-A interaction [30], 
with 184 neutron and 126 proton single-particle wave 
functions and pairing strengths Vo,n = 258.96201 MeV 
and Vo,p = 270.08200 MeV for neutrons and protons, 
respectively. The PES for 240 Pu presents a prominent 
local minimum, corresponding to a fission isomer. If the 
initial test wave functions (harmonic oscillator states in 
the case of SKY3d) are chosen to be prolate with simi¬ 
lar deformations, shape unconstrained static calculations 
converge directly into the isomeric state. This provides 
two initial points for the CHF calculations, starting at 
either the ground or the isomeric state. Some properties 
of the ground state and isomer are presented in Table 
I. These compare well with previous results in the liter¬ 
ature [9]. We consider thus that our choice of pairing 
interaction is suitable. We do not take into account pair¬ 
ing correlations in the dynamical evolution, and hence 
we do not perform a systematic study of the role of the 
pairing interaction. 

One-dimensional projections of the quadrupole- 
constrained PES for 240 Pu are shown in the top three 
panels of Fig. 2. Panel (a) focuses on the dependence 
of the energy on the constrained quadrupole degree of 
freedom. Panels (b) and (c) give the corresponding 
octupole and liexadecapole coordinates associated with 
each quadrupole configuration. For guidance, the color 
scheme is chosen to display changing quadrupoles. Two 
fission barriers are found in the quadrupole degree of free¬ 
dom (top left panel), peaking at Ajo ~ 0.50 and 0.86 re¬ 
spectively. The ground and isomeric states correspond 
to the two minima in the PES next to these barriers. 

Table II presents a comparison between various fea¬ 
tures of the PES calculated in this work and previous 
literature. The table also contains measurements of the 
fission barrier heights and energy differences between 
the ground and the isomeric states, as defined in Fig. 
1. We compare our calculations to recent Hartree-Fock- 
Bogoliubov (HFB) calculations using the Gogny interac- 
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FIG. 2. (Color online) Resulting PES for 240 Pu following a constraint of the quadrupole deformation parameter (? 20 - The top 
panels display the dependence of the constrained energy on the (a) quadrupole, (b) octupole and (c) hexadecapole moments. 
Once the second barrier is overcome, a significant octupole deformation (corresponding to mass asymmetry) develops. Bottom 
panels (d)-(k) display two-dimensional (2D) slices of the 3D density for increasing quadrupole deformation. The isolines 
correspond to 0.05 particles/fm 4 . It is interesting to note that, for states far beyond the second barrier, scission has not yet 
occurred. 


tion [41] , to older HF calculations employing the Skyrme 
Sill parametrization [9] , and to macroscopic-microscopic 
calculations based on the shell-corrected Finite Range 
Liquid-Drop Model (FRLDM) [10]. There is a general 
agreement between the barrier geometries in the theoret¬ 
ical methods. The older Sill calculations predict a large 
second barrier height, most likely because of the assump¬ 
tion of axial symmetry. HF calculations yield barriers 
higher than either the experimental or the macroscopic- 
microscopic predictions, perhaps caused by a lack of dy¬ 
namical effects. The details of the barrier shapes should 
be sufficiently well reproduced in the present calculations 
to deal with the fast fission dynamics beyond the second 
barrier, as well as to give a qualitative description of the 
between-barrier dynamics. 

Panel (b) of Fig. 2 shows a prominent octupole de¬ 
formation setting in at the second fission barrier, as 
would be expected [38]. The relationship between the 
quadrupole and octupole deformation parameters for 
the configurations along the PES is explored in Fig. 3 
(top panel). Although the calculations have been per¬ 
formed constraining only one deformation degree of free¬ 


dom, the observed behaviour is typical for the optimum 
static fission pathway obtained in quadrupole-octupole- 
constrained deformation surfaces calculated using DFT 
[14, 43]. Beyond the second fission barrier, octupole de¬ 
grees of freedom are explored. Panel (b) of Fig. 3 dis¬ 
plays the relationship between the quadrupole and hex¬ 
adecapole deformation parameters. Near the peaks of the 
first and second barriers, the hexadecapole deformation 
sharply drops, and recovers subsequently. This corre¬ 
sponds to a transitioning shape as the neck region of the 
nucleus thins. 

The 3D calculations verify that triaxiality is explored 
at the first fission barrier. The bottom panel of Fig. 3 
shows unambiguously a region of nonzero values of 7 . 
Access to these additional degrees of freedom lowers the 
calculated barrier height with respect to axially symmet¬ 
ric calculations [38]. We note that triaxiality is explored 
significantly in the range 0.36 < /3 2 o < 0.59, but it is 
virtually negligible elsewhere. 

The slices of the density in the lower panels of Fig. 2 
display an increasingly deformed shape as the quadrupole 
degree of freedom grows. Interestingly, the nucleus has 
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E a 

[MeV] 

Eb 

[MeV] 

Eh 

[MeV] 

Method 

Reference 

8.25 

7.68 

3.04 

HF+BCS, Skyrme SkM* 

This work 

8 

13 

4 

HF+BCS, Skyrme Sill 

Table 2 of Ref. [9] 

9.30 

8.40 

3.10 

HFB, Gogny DIM 

Fig. 5 of Ref. [41] 

5.99 

4.91 

2.94 

Shell-Corrected FRLDM 

Table I of Ref. [10] 

6.1T0.3 

6.0T0.50 

2 . 1 T 0.6 

Experiment 

Fig. 27 of Ref. [10] (Madland) 

5.6±0.2 

5.1T0.20 

2.4T0.3 

Experiment 

Fig. 27 of Ref. [10] (Madland) 


TABLE II. Comparison of properties of the fission barrier for 240 Pu from different calculations (defined in Fig. 1). In addition 
to our work, we present the calculations of Flocard et al. [9], Rodrfguez-Guzman and Robledo [41], Moller et al. [10] and the 
experimentally inferred data presented in Ref. [10]. 



FIG. 3. (Color online) (a) Octupole deformation as a func¬ 
tion of quadrupole deformation and (b) hexadecapole defor¬ 
mation as a function of quadrupole deformation for the calcu¬ 
lated PES of 24 () Pu. The octupole deformation, corresponding 
to a mass asymmetry, rapidly onsets after the second fission 
barrier is passed. Panel (c) displays the corresponding 7 de¬ 
formation parameter. 


not fissioned in the range of /3 20 considered, even for 
states beyond the second fission barrier. Despite the 
emergence of a competing fission pathway, a large selec¬ 
tion of increasingly deformed states have been obtained. 
These are useful as starting points for our time-dependent 
calculations. We note that our configurations explore a 
range of states, from configurations with a quadrupole 
deformation less than that of the global HF minimum, 
to configurations well beyond the second fission barrier. 

When performing CHF calculations, we noticed that 
beyond a quadrupole deformation of /3 2 q = 1.25, the con¬ 


figuration jumped abruptly to a competing two-fragment 
fission pathway. This behaviour, owing to the numerics 
as the calculations converged, proved to be unavoidable. 
This can be explained by considering the density slices 
presented in Fig. 4. The masking procedure that is im¬ 
plemented in the calculations around the one-fragment 
configuration does not inhibit a transition to the two- 
fragment configuration, as the two-fragment configura¬ 
tion fits inside the one-fragment masking region. 

The competing two-fragment fission pathway was ex¬ 
plored, starting from the state after the calculations 
jumped pathways. From this configuration, the deforma¬ 
tion was incrementally reduced. This competing path¬ 
way is shown in Fig. 4 with solid squares, and may 
be compared to the original, one-fragment pathway in 
solid circles. Once the quadrupole deformation param¬ 
eter is reduced below /3 20 = 1.01, the HF minimum 
jumps back onto the original fission pathway. The one- 
fragment pathway is also sometimes denoted as fission 
valley, whereas the two-fragment pathway [41, 52] has 
also been referred to as the fusion valley [53, 54]. 

The competing pathway, referred to hereafter as the 
two-fragment pathway [41], displays remarkably differ¬ 
ent configurations to that of the one-fragment pathway. 
Even with identical quadrupole deformations [see panels 
(d) to (g)], the octupole and hexadecapole deformations 
and total energy differ significantly. It is exactly this be¬ 
haviour that the authors of Refs. [10] and [48] identify 
as a flaw when using CHF to explore the PES. Here, we 
exploit this feature to gain an insight on the competing 
fission pathway without having to include a higher num¬ 
ber of constraints in the CHF calculations. 

The fragments in the two-fragment pathway do not 
have an integer particle number. For example, for the 
case of /? 20 = 1.19, the fragments have A 1 = 107.14, 
Z x = 43.14, and A 2 = 132.85, Z 2 = 50.85. We note, 
however, that all the fragments in the two-fragment path¬ 
way correspond, to the nearest integer particle number, 
to 43 7 Tc and si 3 Sb. It would be instructive to project 
the individual fragments onto a good particle number 
[55], because this would give access to a mass distribu¬ 
tion. Our focus here is on the dynamics, though, so we 
postpone this for future work. 
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FIG. 4. (Color online) Panels (a) —(c): One-fragment (solid circles) and two-fragment (solid squares) fission pathways for 
240 Pu. The arrows in panel (a) show the direction in which the PES is explored. Beyond /?20 = 1-25 the one-fragment pathway 
jumps into the two-fragment pathway, and this state is used as the initial configuration for investigating the latter pathway. 
Sample density slices on the competing pathways with the same /I 20 are shown in the panels (d)-(e) and (f)-(g). The isolines 
are separated by 0.05 particles/fm 3 . 


III. TIME EVOLUTION OF CONSTRAINED 
HARTREE-FOCK STATES 

The time evolution of the CHF states obtained for 
240 Pu may be investigated using the TDHF method. 
This analysis will focus on the states on the one-fragment 
fission pathway, starting from configurations beyond the 
fission isomer (that is, those with /3 2 o > 0.68). We in¬ 
vestigate the effect of releasing the imposed shape con¬ 
straints and time evolving the constrained static states, 
to gain an insight of the deformation-induced fission 
(DIF) process in TDHF. 

Our TDHF calculations are performed in a larger grid 
than that used to calculate the initial (static) state. In 
most calculations, the Slater determinant solution to the 
CHF problem is placed in a grid of dimension (42 3 ) 
points, ranging from —20.5 to 20.5 fm in the x, y, and 
2 directions. The BCS occupations associated with each 
wave function remain frozen throughout the dynamical 
simulation. This approximation is relatively harsh, but 
provides a substantial computational advantage. Our 
study represents a preliminary attempt to gauge the po¬ 
tential of TDHF techniques in pre- and post-fission dy¬ 
namics. As such, while time-dependent pairing effects 
are likely to be relevant for quantitative predictions, we 


stay at a relatively qualitative level. Comparisons with 
observables should therefore be considered with care. We 
refer the reader to Refs. [36, 56, 57] for details on more 
thorough treatments of dynamical pairing effects. 

During our calculations, nuclear fragments above par¬ 
ticle emission threshold will be created. Some single par¬ 
ticle wave functions will thus be free to explore the entire 
space of the calculation, up to the box boundaries. Our 
periodic boundary conditions thus may cause artificial 
effects, especially in the analysis of fragment vibrations. 
While spherical TDHF calculations can be performed in 
an analytic continuum [58, 59], the available numerical 
methods that can be applied to 3D cases are compu¬ 
tationally expensive and may not be suitable for large 
amplitude processes [60, 61]. When analysing nuclear 
dynamics, we apply a spatial mask to ensure that ob¬ 
servables correspond only to the nuclei and not dripped 
particles. For separated fragments, individual co-moving 
masks are used. Further details can be found in Ref. [50]. 

As one increases the /3 2 o deformation of the t = 0 state, 
different dynamics take place. Before the second bar¬ 
rier, no fission can occur without collective tunnelling 
not open to the TDHF method. For these initial config¬ 
urations, we see large-amplitude collective motion, which 
has its own interest [50]. This suggests that the TDHF 
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wave functions are exploring a local minimum in multi¬ 
dimensional deformation space and that there is an in¬ 
hibition in rearranging substantially the nuclear density 
while keeping the total energy constant. As a matter of 
fact, the corresponding power spectra are in qualitative 
agreement to those obtained in a giant-resonance calcu¬ 
lation. We do not explore these issues further here, as 
our emphasis is on final fissioning states. 


A. States Between the One- and the 
Two-Fragment Pathways 

Beyond the peak of the static fission barrier, the time 
evolution of several increasingly deformed initial states 
fail to display fission within 9000 fm/c. The evolution of 
the multipole deformation parameters for these states is 
presented in Fig. 5. 

Qualitatively, one observes dramatically different be¬ 
haviour in the time evolution of the multipole deforma¬ 
tions compared to a giant-resonance-like behaviour be¬ 
low the second barrier. For most cases, the elongation 
is seen to increase rapidly during the first 300 to 500 
fm/c, indicated by an increase of /? 2 o- The most ex¬ 
treme case is seen in the bottom left panel of the figure, 
where the quadrupole deformation quickly increases from 
P20 ~ 1-07 to P20 ~ 1-11. 

Beyond the initial increase in elongation, Fig. 5 dis¬ 
plays slow, large amplitude oscillations setting in. Com¬ 
pared to the states below the static fission barrier, these 
oscillations are substantially slower. They therefore 
correspond to lower energy modes in the power spec¬ 
trum. For the initial configurations with /32o = 0.89 and 
0.95 (panels (a) to (f) in Fig. 5), the behaviour of the 
quadrupole deformation is more complex than the other 
two cases (foo = 1-01 and 1.07). The evolution of the 
quadrupole deformation for these cases (fiw = 0.89 and 
0.95) shows a region of rapid increase, then an oscillation 
about a plateau, then another rapid increase followed by 
another plateau. 

An octupole deformation (depicted in the central col¬ 
umn of Fig. 5) is observable in all cases. This is unsurpris¬ 
ing in itself, as the initial configurations are significantly 
octupole deformed. However, an interesting feature is no¬ 
ticeable for the evolution of the states with initial defor¬ 
mation /3 2 o = 0.95 and 1.01. The changes in octupole de¬ 
formation are roughly in phase with either the evolution 
of the hexadecapole parameter or both the quadrupole 
and the hexadecapole parameters. This feature, in addi¬ 
tion to the other differences observed, suggests that the 
mechanism driving the dynamics of the multipole defor¬ 
mations above and below the fission barrier is different. 

The slow, large-amplitude oscillatory behaviour of the 
multipole deformation parameters (second row from bot¬ 
tom of Fig. 5) suggests that, owing to the Coulomb re¬ 
pulsion between the upper and lower lobes, the nucleus 
is attempting to fission. This is in line with the macro¬ 
scopic model of Bohr and Wheeler [62] , where the effect 


of the charge on an incompressible liquid drop is a cru¬ 
cial ingredient to describe the fissioning process. Within 
macroscopic liquid drop models, the surface term com¬ 
petes with the repulsive Coulomb force to inhibit fission; 
it costs energy to form an increasingly deformed shape. 
The TDHF calculations present a similar behaviour, but 
the mechanism arises microscopically. Further calcula¬ 
tions, not shown here for brevity, confirm the relevant 
role of Coulomb repulsion in the dynamics of these con¬ 
figurations [50]. 

Overall, there is significant evidence concerning the 
mechanism responsible for the slow, large-amplitude os¬ 
cillations of the nuclear shape observed on non-fissioning 
configurations beyond the static fission barrier. The 
oscillations are driven by a competition between the 
Coulomb force, trying to cause fission, and the attractive 
nuclear potential terms in the energy functional, counter¬ 
ing this effect. When the states are evolved in time, they 
explore significant collective motion in an attempt to find 
a pathway towards fission. It can only be speculated as 
to whether, with a long enough time evolution, the states 
would eventually fission. Panel (j) of Fig. 5 (evolution of 
a state from /3 2 o=l-07) shows that the quadrupole defor¬ 
mation oscillates around a gradually increasing average. 
This suggests a final state that could eventually fission. 
The timescale for the quadrupole increase is very slow, 
however. A TDHF calculation that explores the time evo¬ 
lution beyond t ss 10, 000 fm/c may begin to encounter 
numerical instabilities. 


B. Intersection of the One and Two-Fragment 
Fission Pathways 

One striking characteristic of the static one- and two- 
fragment pathways is the intersection point that sepa¬ 
rates the initial states which fission upon time evolution 
from those which do not. Figure 6 zooms into the area 
of the PES which is relevant for these differences. We 
divide the PES into three regimes. For configurations 
starting with a deformation below the static fission bar¬ 
rier (/?20 < 0.86), tunnelling is required to reach a fis¬ 
sioned state. This is forbidden in TDHF calculations, in 
which collective coordinates behave semiclassically. We 
therefore define a forbidden region below the barrier. Be¬ 
yond the barrier, in contrast, there is a region where 
fission is inhibited. It will shortly be demonstrated in 
Sec. IIIC that, beyond the intersection of the one and 
two-fragment pathways, fission is allowed within the con¬ 
sidered time scales of the TDHF calculations. The line 
between /3 2 o = 1-07 and 1.10 thus separates the inhibited 
and allowed regions for fission, but the deformation itself 
does not correspond to a threshold in the dynamic calcu¬ 
lations. The bottom row of Fig. 5 indicates, for example, 
that the state with an initial deformation just below the 
separating line can evolve dynamically to a state with 
deformations beyond this very same line, but without 
fissioning. 
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FIG. 5. (Color online) Time evolution of quadrupole (left columm), octupole (centre columm) and hexadecapole (right columm) 
multipole parameters for initial states which are solutions to CHF calculations, with initial quadrupole deformation labelled 
on the right-hand side. All of the initial states are deformed beyond the peak of the static fission barrier, P 20 > 0.86. 


An intuitive explanation may be given for the signif¬ 
icance of the intersection point of the pathways on the 
PES with regard to fission occurring upon time evolution. 
In TDHF, the total energy is conserved. It is because of 
the inclusion of excitations (internal and translational) 
that nuclear configurations may change upon time evo¬ 
lution. For the states which undergo fission (allowed re¬ 
gion in Fig. 6, 0 2 0 > 1.085), for a given value of 0 2O , the 
two-fragment state is more bound than the one-fragment 
state. In consequence, the one-fragnrent state can evolve 
into a two-fragment configuration at a constant 0 2 q by 
releasing nuclear potential energy into excitation energy. 
Of course, the picture is not really that simple, because 
the significance of the static PES becomes less clear in 
the dynamic case. In general, configurations which do 
not correspond to the static fission pathways will be ex¬ 
plored. Further, a slight change of configuration will be 
required to move from the static one-fragment state to a 
fissioned configuration. In other words, the exact config¬ 
urations on the two-fragment pathway cannot be reached 
dynamically from the one-fragment pathway, but an ex¬ 
cited two-fragment configuration of a similar deforma¬ 
tion can. The reasoning presented is that, as the static 
two-fragment state is more bound than the correspond¬ 
ing one-fragment state, the optimum TDHF trajectory 
is to evolve the one-fragment static state towards an ex¬ 


cited fissioned configuration by undergoing only a modest 
rearrangement of the nuclear shape. 

In contrast, in the inhibited region of Fig. 6 (0.86 < 
P20 < 1.085), for a given 0 2O in the one-fragment path¬ 
way, the two-fragment state with the same 0 2 q is less 
bound. Owing to energy conservation, the one-fragment 
state cannot move to the two-fragment state at the same 
0 2 q. The only way to reach a two-fragment solution of 
equal binding energy (or an excited configuration with 
greater binding energy) is through a significant change 
in deformation and rearrangement of the nuclear state, 
which accounts for the inhibiting time scale for fission 
to occur. It may be that the presented calculations 
lack the degrees of freedom necessary to allow a fission¬ 
ing path to be found in this window at all, and that a 
method beyond basic TDHF, including either collisions 
[63] or dynamical pairing effects [36], is needed to reach 
the fissioned configuration. Further investigation of the 
link between static and dynamic configurations apply¬ 
ing density-constrained TDHF [35, 64], would certainly 
be of interest. This method allows the dynamic config¬ 
urations to be “frozen”, removing internal and collective 
excitations, thus bridging between static and dynamic 
configurations. 

In our frozen BCS approach, however, the lack of dy¬ 
namical single-particle occupation effects could be par- 
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FIG. 6. (Color online) Static one- and two-fragment fission 
pathways (solid circles and squares, respectively). For TDHF 
calculations, three regions may be defined when following the 
one-fragment pathway. In a first region, fission is forbidden 
within TDHF time scales. Beyond the maximum in the PES 
at P 20 = 0.86, the time scale for fission is inhibiting. Finally, 
beyond the crossing of two pathways at P 20 = 1.085, fission 
is allowed. See text for more details. 


1 

Forbidden 

1 1 

Inhibited 

n 1 

Allowed 

" 

■ 




_ 

. 


1 1 1 

f 

• 

f 

f ^ 

T 

__ 

_1_ 

IF • 

2F m 

_ 1 _ 1 _ 

_ 

j_ J _ 


.7 0.8 0.9 1 1.1 1.2 1 

P20 


ticularly important in this region of the PES. Relatively 
compact initial states can have substantially different 
single-particle structure as compared to the real fission¬ 
ing fragments. The dynamical final state that we gener¬ 
ate is therefore only an approximation to the real one, 
which will be captured better in simulations including 
time-dependent superfluidity [36]. While the dynamics 
in the interface between the inhibited and allowed re¬ 
gions can be considered somewhat artificial, we expect 
these effects to be less important for initial states well 
into the allowed regime, foo > 1.15. There, the struc¬ 
ture of the pre-fission system is already reminiscent of 
a two-fragment state, and dynamical rearrangement of 
single-particle orbits is likely to be less important. 


C. Fissioning States 

For the static states with a quadrupole deformation at 
and beyond the threshold of /3 2 o = 1 . 10 , binary fission is 
seen to occur as the wave functions are evolved in time. 
The calculations to obtain the data for this Section were 
performed in a larger grid of size 48 x 48 x 160 points, 
corresponding to —79.5 to 79.5 fm in the 2 direction, 
and —23.5 to 23.5 fm in the x and y directions. The 
calculations were set to end once the separation of the 
center of mass of the two fragments exceeded 100 fm. 
This cutoff avoids spurious effects associated with the 
fragments approaching the grid boundaries. 

Figure 7 shows the typical time evolution of the parti¬ 
cle density for the fissioning case by presenting 2D slices 
of the 3D density at various times for the state with 
initial deformation /3 2 q = 1.19. It is difficult to estab¬ 


lish exactly the scission point in a calculation involving 
quantum mechanical wave functions and densities. We 
take an operational approach and define “scission” as 
the point when the minimum density between the frag¬ 
ments along the principal axis of the system is less than 
0.05 particles/fm 3 . As we shall see in the following, this 
also corresponds to the point where a sizable collective 
energy develops as the fission products begin spatially 
separating. For the case where the initial quadrupole de¬ 
formation is P 20 = 1-19 (presented in Fig. 7), it takes 
between 775 and 800 fm/c for the density between the 
two fragments to drop below this threshold. Figure 8 
displays sample current vectors corresponding to the par¬ 
ticle density slices of Fig. 7. The current vectors display 
the system smoothly transitioning into a two-fragment 
configuration. There is no dramatic rearrangement of 
the density during time evolution. Throughout the cal¬ 
culation, the currents in the two preformed fragments 
are clearly distinguishable, and do not interact with one 
another. The central region has negligible current, and 
the two lobes stretch against each other. The magni¬ 
tude of the current vectors in Fig. 8 gradually increases 
as the fission occurs. Beyond the point of scission, they 
increase rapidly as the fragments accelerate away from 
one another. 

States with increasingly large initial deformations fis¬ 
sion in a qualitatively similar way. We now proceed 
to compare their fission outcomes in terms of macro¬ 
scopic observables. The states with static deformation 
/3 20 = 1.10,1.13,1.19 and 1.25 were evolved in time to 
investigate the fission of the different initial configura¬ 
tions. The time evolution of the multipole moments for 
these states are shown in Fig. 9. Different nuclear shapes 
are explored as the nucleus evolves from the various static 
states. Other than the case with initial /3 2 0 = 1.25, we 
find that as /? 2 o and ,840 increase, the octupole defor¬ 
mation, as reflected in ^ 30 , remains virtually constant. 
Here and in the following discussion, we sharply cut off 
the time evolution at the point of scission. An analysis 
of the post-scission fragments is presented in Sec. IV. 

The chosen dynamic pathway towards fission, depend¬ 
ing on how the particles rearrange during the time evolu¬ 
tion, may have significant consequences upon the proper¬ 
ties of the post-fission system. This will produce a range 
of fission fragments dependent on the initial configuration 
which is time evolved. Once again, this differs from the 
static case, where CHF calculations following the one- 
fragment fission pathway will only produce one resulting 
fissioned configuration. The distribution of fission prod¬ 
ucts obtained with TDHF is in line with experimental 
investigations (see Sec. IV). 

The time scale required for the initial configuration to 
fission varies. The least elongated case, with /3 20 = 1.10 
takes tscission ~ 1250 fm/c for scission to occur. Figure 9 
shows the quadrupole deformation increasing rapidly for 
approximately 600 fm/c [panel (a), solid line]. Between 
600 and 1200 fm/c, the rate of increase in quadrupole 
deformation reduces as the nucleons rearrange out of the 
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FIG. 7. (Color online) Slices of the total particle density for various times, starting from the static case with (3 20 = 1.19. It 
takes between 775 and 800 fm/c for scission (as defined in the text) to occur. The isolines are separated by 0.05 particles/fm 3 . 


neck into the upper and lower fragments. Small oscilla¬ 
tions in the octupole deformation can be seen as the sys¬ 
tem transitions into the preferred configuration. Beyond 
1200 fm/c, the neck rapidly vanishes as the fragments 
take form and begin to separate, resulting in an accelera¬ 
tion in the increase of the quadrupole and hexadecapole 
parameters. 

For more deformed initial states, the time taken to fis¬ 
sion is significantly shorter. This can be explained as the 
initial configuration has fewer particles in the neck region. 
Upon time evolution, less rearrangement is required for 
the two fragments to take form, and the Coulomb interac¬ 
tion rapidly drives the configuration to fission. The most 
extreme case investigated here is that with initial defor¬ 
mation P 20 = 1.25. Panel (k) of Fig. 2 shows the initial 
density for this state. Two fragments are already taking 
form, connected only by a thin elongated neck. This nar¬ 
row structure rapidly dissipates into the top and bottom 
fragments upon time evolution, and within t sc i ss ion ~ 400 
fm/c the system fissions. 

The evolution of the decomposed contributions to the 
energy density functional for the fissioning cases are pre¬ 
sented in Fig. 10. The decomposed energy density func¬ 
tional for the entire two-fragment system is shown. The 
total energy of the two separate fission fragments will 
be analysed separately in Sec. IV. Figure 10 displays the 
evolution of the energy density functional up to and be¬ 
yond the point of scission. The total energy, shown on 
the panel (i), is, in principle, conserved within the TDHF 
calculations. For these fission calculations, the fluctua¬ 


tions during the time evolution are less than 4 MeV. 

The dynamic calculations allow for both translational 
motion and internal excitations. In the fissioning case, 
the nuclear binding energy is expected to be transformed 
mainly into the translational kinetic energy of the frag¬ 
ments. The nuclear collective kinetic energy is conven¬ 
tionally defined within TDHF as, 

h 2 f j(r ) 2 

-Ucoll. kin. = X- / dr , (2) 

2 m J p(r) 

where p(r) is the particle density and j ( r) is the current 
density [30]. This collective kinetic energy contains con¬ 
tributions from internal excitations of the nucleus, such 
as resonances, as well as the translational kinetic energy 
of the post-fission fragments. It is presented in panel (h) 
of Fig. 10. The total kinetic energy is shown separately, 
in panel (g). It is difficult to untangle the collective exci¬ 
tation energy attributed to the internal excitation of the 
fission fragments, from that attributed to translational 
motion. In fission reactions, the energy release is typi¬ 
cally attributed to be « 80 % in the form of translational 
energy, and the other « 20 % is released in the form of 7 
rays, prompt neutron emission, and radioactive decays of 
the fragments [65]. In our TDHF calculations, only part 
of these effects can be described. The degrees of freedom 
to allow fission fragments in hot resonance states to de¬ 
cay by particle emission, for instance, are included. We 
demonstrate that the excitation energy of the fissioned 
system in the TDHF simulations is dominated by the 
























11 


t=0 fm/c 


t=250 fm/c 


t=500 fm/c 


t=600 fm/c 


15 

10 

I 0 

N .5 
-10 
-15 


15 

10 

I 0 

N -5 
-10 
-15 


t=650 fm/c 


“i—i—i—i—i—i—r 

IT 

VVtV 
! fVj 


t=750 fm/c 



t=675 fm/c 




n—i—i—i—i—i—r 


wiv ;* 1 


l-M 1 


i__i_i_L_ 


t=850 fm/c 


1—i—i—i—i—i—r 



x [fm] 


x [fm] 


x [fm] 


x [fm] 


FIG. 8. (Color online) Current vectors corresponding to the slices of the particle density presented in Fig. 7. The vectors have 
been normalised to the same scale in each panel, to a visually instructive length. Note that the dimensions of the grid exceed 
what is presented in these panels. 


translational kinetic energy, with a small contribution 
from internal collective excitation of the fragments. This 
will be discussed in Sec. IV C. 

During the time evolution, the individual components 
of the energy functional may be examined separately. 
The physical interpretation of the evolution of each term 
of the integrated energy functional shown in Fig. 10 may 
not necessarily be simple. It is useful to identify which 
densities contribute to the separate terms to qualitatively 
explain the behaviour of the energy functional: 

• Eq and E 3 terms: These terms, proportional to the 
Skyrme parameters to and / 3, are the central terms 
of the functional. They are attractive and repul¬ 
sive,! respectively [see panels (a) and (d)], pro¬ 
viding a similar functional form with an expected 
cancelation between them. They depend upon the 
particle density. As the nucleus approaches scis¬ 
sion, the Eq term is reduced in strength, thanks to 
the small, unfavourable density that briefly exists 
in the neck region. There is then a sudden increase 
in biding as separate, more stable fragments are 
formed. After this, the energy contributions os¬ 
cillate with a smaller magnitude than the changes 
during the fission process, corresponding to the ex¬ 
cited collective motion of the fragments 


• Ei term: The E\, shown in panel (b), term con¬ 
tains contributions from the kinetic, particle, and 
current densities. The time evolution of this term 
is qualitatively similar to the Eq and E 3 contribu¬ 
tions, with the same sign as the E 3 term but an 
absolute smaller scale. This suggests the that den¬ 
sity that governs the E 0 and E 3 terms, the particle 
density, is also the most relevant contribution driv¬ 
ing the Ei term. 

• E 2 term: The E 2 contribution in panel (c) contains 
the Laplacian of the particle density, and is com¬ 
monly associated with a surface term. As the par¬ 
ticles rearrange into the two fission fragments, this 
term increases in magnitude - i.e. gives an overall 
more repulsive contribution to the entire system. 
This can be explained by the two-fragment system 
having a combined surface fraction which is greater 
than that of the initial configuration. The gain in 
energy for this term up to the point of scission is 
dependent upon the deformation of the initial con¬ 
figuration. It can increase by as much as 45 MeV 
for the static configuration with P 20 = 1-10. 

• Coulomb term : The Coulomb energy depicted in 
panel (f) is determined from the distribution of the 
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FIG. 9. (Color online) Time evolution of (a) quadrupole, 
(b) octupole and (c) hexadecapole deformation parameters 
for various static configurations observed to fission upon time 
evolution. The evolution is stopped at the scission point, as 
defined in the text. 


charged protons. The magnitude of the repulsive 
Coulomb term slowly decreases as the nucleus elon¬ 
gates. The overall reduction is of the order of 200 
MeV as the system evolves. At the point of scis¬ 
sion, the rate at which the term reduces rapidly 
accelerates, as the two charged fragments separate 
from one another in co-ordinate space. At infinite 
fragment separation, the Coulomb term will reduce 
to the contributions of the Coulomb energy for each 
nucleus, without further interactions. 

• Kinetic and Collective Kinetic terms : The kinetic 
energy can be determined from integrating the ki¬ 
netic density. As mentioned above, the contribu¬ 
tion to this energy from collective motion (assumed 
to be predominantly translational beyond scission, 
rather than internal collective excitation) can be 
decomposed according to Eq. (2) and we show it in 
panel (h). The collective energy is initially small, 
as it is only associated with the internal currents 
as the nucleus slowly rearranges into a fissioned 
configuration. The state with initial deformation 
P 20 = 1-10 shows the most gradual transition to 
fission. An initial rapid increase and a satura¬ 
tion in the collective energy is seen before scission, 
which corresponds to the (previously discussed) 
rapid initial elongation, then extended rearrange¬ 
ment phase as the configuration evolves (see Fig. 
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FIG. 10. (Color online) Evolution of the decomposed energy 
density functional for the fissioning systems. The calculations 
are terminated when the fragments are separated by 100 fm. 
For reference, the vertical lines in the panels corresponding to 
kinetic energies show the corresponding scission points. The 
total energy is conserved within fluctuations no greater than 
4 MeV. 


9). In contrast, the state with /? 2 o = 1-25 is already 
close to the point of scission, so that the Coulomb 
interaction between the two lobes rapidly drives 
the configuration to the scission point (within the 
first few hundred fm/c), where translational motion 
rapidly accelerates once the neck ruptures. This 
shorter timescale could explain the more extreme 
behaviour observed in the evolution of the other 
terms in the energy functional as the particles in 
the neck have less time to rearrange into the two 
fragments. At the point of scission, the collective 
kinetic energy rapidly increases at a rate similar to 
that of the reduction in the Coulomb energy. The 
threshold collective kinetic energy associated with 
the scission point is between 6 and 8 MeV in all the 
cases presented. The vertical lines corresponding to 
the scission points are displayed in panels (g) and 
(h) of Fig. 10. The definition adopted for scission 
[66] is justified by considering the rapid increase of 
collective kinetic energy at this point. The gain in 
the total kinetic energy beyond the scission point 
can be attributed to the gain in collective energy, 
and is of the order of 150 MeV in the time consid¬ 
ered. This relates to the loss in Coulomb energy 
beyond the point of scission, as would be expected. 
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• Spin-Orbit term: The spin-orbit potential is im¬ 
portant for the single-particle structure. Spin-orbit 
partner levels combine, when fully occupied, to give 
a total zero energy contribution, so the contribu¬ 
tions reflect the changing of the single particle lev¬ 
els and the change of meaning of spin-orbit partners 
in the parent nucleus to the daughter nuclei. One 
therefore expects this contribution to be relevant in 
determining the structural details of the final fis¬ 
sion fragments. The final approximately constant 
values observed in panel (e) following scission cor¬ 
respond to the sum of the two independent spin- 
orbit terms of the separate fragments. The notably 
different behaviour of the term for the state with 
initial deformation /3 2 o = 1-25, compared to the 
others, suggests that different shell effects are act¬ 
ing. Indeed, the masses of the fission products are 
significantly different for this case (see Sec. IV A). 
The spin-orbit term also has significant contribu¬ 
tions from time-odd densities, so that the evolution 
can explore configurations which may not be acces¬ 
sible on the static PES. Overall, the term varies by 
less than 30 MeV during time evolution. 


In experimental studies of fission, it is customary to 
measure the kinetic energy of the fission fragments. We 
can find an analogous observable within TDHF by mak¬ 
ing use of the collective kinetic energy, defined by Eq. 
(2), and assuming that the translational kinetic energy 
dominates this term. Referring to Fig. 10, the collective 
kinetic energy and Coulomb energy are both expected to 
plateau as the separation of the two fragments becomes 
large. Unfortunately, the Coulomb force is long-ranged. 
Because increasing the dimensions of the numerical grid 
is extremely computationally expensive, we take an alter¬ 
native perspective. We extrapolate the collective kinetic 
energies to large times to estimate the asymptotic value 
of the collective energy as the separation r tends to oo. 

A simple approximation to the time-dependent be¬ 
haviour can be made from classical mechanics. Let us 
assume two pointlike fragments with charges Z u and Z\. 
and masses M u and M/. If the two fragments fission 
from a ground state because of the Coulomb force, and 
convert all this energy into translational kinetic energy, 
energy conservation implies: 

\M u vl + \Miv1 = • (3) 

Here, Mj, Vi and Zi are the mass, velocity and charge 
of each fragment (i = u for upper and i = l for lower 
fragment). The constant n is the Coulomb constant. As 
momentum must be conserved, 


M u v u + Mivi= 0, 

Eq. (3) may be rewritten, substituting for v v 

Z„Z, 


Vl 


M? 

jt, +Ml]=2Kl 


(4) 

(5) 


For a given fissioned system, M u , Mi, Z u and Zi are con¬ 
stant. A differential equation for ^ (= vf) can be formed 



where all the constants are combined into 0. Performing 
the integration 



allows the solution 

r 3 / 2 =r- 3 0 /2 + ^Ve(t-t 0 ) (8) 

to be written. According to this approximation, the dis¬ 
tance between the two fission fragments, r, is approxi¬ 
mately proportional to t 2 ^ 3 . By assuming that the loss 
in Coulomb energy is equal to the gain in collective ki¬ 
netic energy (that is , Fcoui = E co n. kin.), a fit of the 
form 

/W = ° + ( t _ 6 c) 3/2 ( 0 ) 

can be performed to interpolate the collective kinetic en¬ 
ergy to asymptotically large values of t. Figure 11 shows 
a sample interpolation of the collective kinetic energy as¬ 
suming the above form for the case of initial deformation 
P 20 = 1.10. The fit is performed over three time ranges: 
once the centers of mass are separated beyond 30, 50 and 
60 fm, respectively. As the separation tends to oo, the 
fit parameter a can be interpreted as the final collective 
kinetic energy. Table III contains the values obtained for 
each of the fissioning cases with different distance fits. 
As a crude method to represent the uncertainty in the fi¬ 
nal collective energy value, we present the mean and the 
standard deviation of the three interpolation fits. 

The values shown in Table III demonstrate that the 
resulting collective kinetic energy varies by about 10% 
depending upon the region of the data the fit is per¬ 
formed to. In this very crude model, fragment defor¬ 
mation, particle emission (discussed in the next section) 
or tidal effects associated with the extended nature of 
the nuclei are not accounted for. This suggests that the 
results obtained from the interpolation method should 
serve only as illustrative values. We note that an alter¬ 
native method, based on computing the collective kinetic 
energy from the center-of-mass momenta of each frag¬ 
ment, provides quantitatively similar results. 

In comparison to the experimentally measured kinetic 
energy of the fissioning systems displayed in Table IV, the 
theoretical values presented in Table III are of a similar 
order of magnitude. The crude theoretical estimates are 
about 5 to 25 MeV larger in all cases. We note, however, 
that the experimental values correspond to an average 
kinetic energy. We only have access to a single fissioning 
event per static state, and a larger sample of theoretical 
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FIG. 11. (Color online) Fits to the obtained collective ki¬ 
netic energy for the initial state with deformation 020 = 1.10. 
Fits are performed over three different ranges: from the point 
where the separation of the fragments exceeds 30, 50 and 60 
fm, respectively. 


TABLE III. Interpolated total kinetic energy corresponding 
to different initial configurations. The fit of Eq. (9) has been 
performed once the fragment separation exceeds 30 fm, 50 fm 
and 60 fm. 


Static 

Deformation 

020 

Coll. KE 
(30 fm fit ) 
[MeV] 

Coll. KE 
(50 fm fit ) 
[MeV] 

Coll. KE 
(60 fm fit ) 
[MeV] 

Mean 

± St. Dev. 
[MeV] 

1.10 

210.3 

206.5 

203.4 

206(4) 

1.13 

210.8 

200.0 

193.7 

202 (8) 

1.19 

205.8 

196.8 

191.3 

198(8) 

1.25 

193.4 

180.8 

176.3 

183(9) 


results would be required to enable a quantitative com¬ 
parison. Further discussion of methods to deduce the 
energy released by the fission reaction within TDHF is 
presented in Sec. IV B. 

TABLE IV. Measured total kinetic energies from experiments 
on spontaneous fission, thermal neutron-induced fission and 
photo-fission. The measurements correspond to the pre¬ 
neutron emission fragment energies. 


Method 

Kinetic energy [MeV] 

Reference 

240 Pu(s.f) 

178.85i0.30 

[2] 

240 Pu(s.f) 

179.00i0.08 

[3] 

239 Pu(rW) 

177.69 

[2] 

239 Pu(n th ,/) 

177.65i0.01 

[3] 

240 Pu( 7 ,/) (12 MeV) 

176.39i0.24 

[2] 

240 Pu( 7 ,/) (15 MeV) 

175.80i0.24 

[2] 

240 Pu( 7 , /) (20 MeV) 

175.15i0.24 

[2] 

240 Pu( 7 ,/) (30 MeV) 

174.98i0.31 

[2] 


IV. FRAGMENT ANALYSIS 

Beyond the point of scission, we consider a two- 
fragment system. The published distribution of SKy3d 
has some capacity to analyse two-fragment dynamics [30] , 
and a version has been modified further to investigate 
the fissioning system and extract some useful observables 

[50]. 

A. Mass Distributions 

As the post-fission fragments are excited, they may 
decay by particle emission. TDHF displays this decay 
by the spreading of the single particle wave functions 
from the region of central density, corresponding to the 
nucleus. When masking the region around the nucleus, 
this decay results in a reduction in the integrated particle 
density over time. For the cases of DIF examined in this 
paper, this decay is of the order of 0.1 — 0.2 particles 
during the whole postscission time evolution. 

To compare to experimental studies, we identify the 
number of particles in each fragment prior to any parti¬ 
cle emission. This is done by integrating the total density 
in each half of the numerical grid separated by the divid¬ 
ing plane immediately after the scission point. The inte¬ 
gral is performed without any masking. An uncertainty 
in the particle number of the fragments may be associ¬ 
ated with the fluctuation of this measurement through¬ 
out time evolution, which is less than 0.05 particles for 
all the considered cases. We are thus confident that the 
fission fragments that we produce have a good average 
mass number. 

These fragment masses can be compared directly to 
experimental data. Table V displays the resulting fission 
fragment masses obtained from this theoretical study. 
We note that these can hardly be considered distribu¬ 
tions, but rather should be taken as indications of what 
dynamical calculations can contribute to fission stud¬ 
ies. The two-fragment static configuration is included 
for comparison. It is important to stress that this differs 
from the deformation-induced fragments. It produces 
more symmetric states, which brings it closer to experi¬ 
mental results of neutron-induced fission. The Table also 
includes the particle number rounded to the nearest in¬ 
teger. It would be of interest to project the individual 
fragments onto the particle number [55] , to obtain a mass 
number distribution, although this would not necessarily 
agree with experiments [36]. 

Further, Table VI contains experimental data taken 
from Refs. [2, 3], listing the most likely masses of the 
fission fragments. The references study various fission 
processes in 240 Pu, including spontaneous fission, ther¬ 
mal neutron-induced fission, and various energy photon- 
induced fission reactions. The spontaneous fission data 
has been included for completeness. 

We emphasise that fission produces a range of masses. 
The values quoted in Table VI correspond only to the 
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TABLE V. Fission fragments obtained from evolving initial static configurations from the one-fragment fission pathway. The 
uncertainties in the particle numbers are a conservative estimate related to the fluctuation in the particle number in the region of 
the grid corresponding to the separate fragments throughout time evolution. The result from the static two-fragment pathway 
for @20 = 1.19 is included for comparison. 


Static 

Deformation /I 20 

Heavy Fragment 
A,Z 

Light Fragment 
A,Z 

Heavy Frag. 
(Integer) 

Light Frag. 
(Integer) 

1.10 

1.13 

1.19 

1.25 

1.19 (2f) 

136.33(5) , 52.78(5) 
135.02(5) , 52.23(5) 
136.13(5) , 52.70(5) 
143.70(5) , 55.65(5) 
132.81 , 50.84 

103.67(5) , 41.23(5) 
104.98(5) , 41.77(5) 
103.87(5) , 41.30(5) 
96.30(5) , 38.35(5) 
107.04 , 43.13 

136t 

53 1 

If Te 

136t 

53 1 

ifCs 

sfSb 

!? 4 Nb 

i!fMo 

ii 4 Nb 

96 c 

38 

i°7 T r 

43 lc 


TABLE VI. Experimentally measured average masses following the fission 24() Pu. The measurements for neutron-induced 
fission were taken before neutron emission of the fissioned fragments. 


Method 

Heavy Fragment 

Light Fragment 

Reference 

240 Pu(s.f) 

138.74i0.20 

101.26i0.20 

[2] 

240 Pu(s.f) 

138.96i0.04 

101.31i0.04 

[3] 

239 Pu(rW) 

139.67 

100.33 

[2] 

239 Pu (n th ,f) 

139.73i0.01 

100.27i0.01 

[3] 

240 Pu( 7 ,/) (12 MeV) 

139.88i0.14 

100.12i0.14 

[2] 

240 Pu( 7 ,/) (15 MeV) 

139.92i0.09 

100.08i0.09 

[2] 

240 Pu( 7 , /) (20 MeV) 

139.84i0.08 

100.16i0.08 

[2] 

240 Pu( 7 , /) (30 MeV) 

139.71i0.14 

100.29i0.14 

[2] 


most likely fissioned configuration. This is a crude com¬ 
parison of experimental data to theoretical calculations. 
A more meaningful comparison would go through a char¬ 
acterization of the full mass distribution, in line with the 
recently proposed method of Ref. [39] . While we cannot 
at present produce a full mass distribution from projec¬ 
tion into particle numbers, we can use each of the four 
different fission cases to build a schematic fission frag¬ 
ment mass distribution. Referring to Fig. 12, which dis¬ 
plays data for neutron-induced fission at low energies in 
panel (a), the obtained theoretical values in panel (b) 
fall well within the experimentally obtained mass distri¬ 
bution. With the limited data set available, the TDHF 
results seem to be consistent with the experimental data. 
Figure 12 also includes for comparison the adiabatic, con¬ 
strained static result of the two-fission pathway. All in 
all, dynamical effects seem to induce a larger mass asym¬ 
metry in the fission fragments. 


B. Energy of Fission Fragments 

By applying masks around the spatial regions of the 
fission fragments, the energy density functional corre¬ 
sponding to each individual fragment may be obtained. 
We note, however, that the interpretation of the results 
is not as simple as in the two-fragment case. The nuclear 
part of the energy density functional is short-ranged. We 
therefore expect that an integral in the spatial region cor¬ 


responding to the individual fragments will be a faithful 
representation of the nuclear part of the energy density. 
The Coulomb interaction, however, is long-ranged. As 
well as the Coulomb interaction within the individual 
fragments, there is a contribution from the interaction 
with one another. This contribution is missing in the 
present two-fragment integral. Further, the fragments 
decay by particle emission, which imparts some time de¬ 
pendence upon the integrated energy corresponding to 
the individual fragments. 

The time evolution of the total integrated energy func¬ 
tional corresponding to the heavy and light (lower) fis¬ 
sion fragments is shown in panels (a) and (b) of Fig. 
13, respectively. The time measurement at t sep = 0 be¬ 
gins when the fragments are sufficiently separated such 
that the masks no longer overlap. To a good approxima¬ 
tion, the fragment total energies are constant over time. 
A slight drift is observed in the time evolution, owing 
to the long-range effects of the Coulomb interaction as 
well as particle decay. We only provide the total inte¬ 
grated energy of the post-fission fragments in Fig. 13. 
The evolution of the decomposed terms corresponding to 
the individual fragments shows no remarkable behaviour. 
Hereafter the fragment energy at the cutoff time is de¬ 
noted as E*. 

The total energies of the two fragments provide an al¬ 
ternative way of computing their collective kinetic ener¬ 
gies. The total excitation energies of the fragments is the 
sum of their translational and internal collective kinetic 
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FIG. 12. (Color online) (a) Experimental independent fission 
yields for neutron-induced fission at E = 0.0253 eV. The data 
are from Ref. [67]. (b) Theoretical mass fragment distribu¬ 
tions. The red bars correspond to the binned TDHF results 
and the blue bar corresponds to the static two-fragment mass 
split. See text for more details. 


energies. If we subtract the total integrated energy to the 
corresponding ground state energy, we obtain a new es¬ 
timate for the fragment excitation energy. This method 
complements the approach of interpolating the total col¬ 
lective kinetic energy of the system, as presented in Sec. 
Ill C. As long as the ground state and the fragment ener¬ 
gies are qualitatively correct, this second method should 
produce comparable results. 

The solver SKy3d has thus been applied to deduce the 
ground states of the fission fragments to the nearest in¬ 
teger particle number. The result for the ground-state 
energy of these isotopes is presented in column 4 of Ta¬ 
ble VII. Here, it is debatable whether the energy func¬ 
tional in SKy3d contains all the terms required to cal¬ 
culate odd-odd and odd-even nuclei. The full time-odd 
contribution is presented in Ref. [68], and the functional 
in SKy3d does not include all these terms. However, as 
the functional used for the static calculations is consis¬ 
tent with that applied to dynamic calculations, Galilean 
invariance is conserved. The functional used in SKY3d 
therefore satisfies all the invariance properties required 
to perform static calculations of odd-odd and odd-even 
nuclei, even if the functional is not in its most “complete” 
form. We also note that pairing has been neglected in 
these calculations. The total pairing contribution to the 
energy functional is typically of the order of 0 — 10 MeV, 
which is small compared to the excitation energies in the 
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FIG. 13. (Color online) Summed energy density functionals 
for the region of space corresponding to the (a) heavy and 
(b) light fission fragments. The quadrupole deformation of 
the initial state is labelled. The drift in the energy can be 
attributed mainly to the Coulomb interaction between the 
two fragments. See text for more details. 


fissioning case. 



Initial P 2 o 


FIG. 14. (Color online) Comparison of the summed A E from 
Table VII (solid circles) to the mean interpolated collective 
kinetic energy from Table III) (solid triangles) for each of the 
fissioning cases. Indicative error bars are displayed. See text 
for details. 


We have access to two measurements of the total exci¬ 
tation energy of the systems, either by interpolating the 
evolution of the collective kinetic energy or by comparing 
the ground-state fragment energies to the excited frag¬ 
ment energies. Figure 14 displays the mean interpolated 
collective kinetic energies (circles) presented in Table III. 
We compare these to the total fragment excitation en- 
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TABLE VII. Comparison of the fission fragment energies to the ground-state energy calculated using the SkM* interaction. 
The fragment total energy at the cutoff time is denoted by E* (see Fig. 13), the ground state energy by E gs , and the difference 
( E * — E ga ) by A E. See text for more details. 


Static o 

Heavy 

Frag. 

E* 

[MeV] 

Egs 

[MeV] 

A E 
[MeV] 

Light 

Frag. 

E* 

[MeV] 

Egs 

[MeV] 

A E 
[MeV] 

1.10 

1.13 

1.19 

1.25 

136t 

53 1 

sfTe 

136t 

53 1 

M 4 Cs 

-1029.22 

-1023.81 

-1034.23 

-1090.17 

-1118.31 

-1110.24 

-1118.31 

-1162.47 

89.09 

86.43 

84.04 

72.30 

ifNb 

i2 5 Mo 

i? 4 Nb 

96 q 
38^ r 

-747.86 

-757.38 

-749.41 

-697.78 

-854.54 

-865.71 

-854.54 

-796.94 

106.68 

108.33 

105.13 

99.16 


ergy (AU hea v y f rag . + AUnght f rag .) for each fissioning case 
(triangles). The error bars in the values of A E display 
an uncertainty of 10 MeV, a conservative estimate for 
typical pairing correlations, possible deformation ener¬ 
gies, and also that only nearest-integer nuclei are consid¬ 
ered. For the interpolated collective energy, the values 
presented are the mean of the three interpolations per¬ 
formed at different fragment separations. We present a 
preliminary, rough estimate of the error of this calcula¬ 
tion, obtained from the standard deviation of these three 
data points (see Table III). 


Figure 14 shows that within the error bars, the results 
from the two techniques produce consistent values of the 
energy released in the fission process. A smaller error 
for both predictions could be achieved by performing the 
calculations up to the point where the Coulomb inter¬ 
action is negligible, and obtaining ground-state energies 
incorporating pairing correlations. As mentioned, the en¬ 
ergies presented here are measures of the total excitation 
energy of the system. We note that this includes both 
excitation energy of the fragments and their relative mo¬ 
tion, in contrast to the experimental definition that often 
neglects the latter. Section IV C will demonstrate a tech¬ 
nique which may be used to decouple the translational 
kinetic energy from the internal collective excitation en¬ 
ergy of the two separate fragments. These will provide 
further insight into the excitation mechanisms in TDHF 
fragment formation. 


Figure 15 shows experimental measurements of the 
kinetic energy reproduced from Ref. [3] for thermal 
neutron-induced fission in 240 Pu. The range of collec¬ 
tive kinetic energies deduced in this section are marked 
with a shaded box. By attributing the deduced total 
excitation energies solely to translational kinetic energy, 
this assumes that the internal collective excitation of the 
fragments are comparatively small. This will be demon¬ 
strated shortly. Despite the limited sample of theoretical 
data, the results agree with the experimental range of 
values. In particular, our results fall well within the ex¬ 
perimental distribution. 



Kinetic Energy [MeV] 


FIG. 15. (Color online) Experimental kinetic energy distri¬ 
bution taken from Ref. [3]. The shaded box corresponds to 
the extremes of the range of kinetic energy values displayed 
in Fig. 14. 


C. Collective Excitation Modes of Fission 
Fragments 

As mentioned, the excitation energy of the fission frag¬ 
ments is assumed to be dominated by translational ki¬ 
netic energy. However, as well as translational motion, 
the fragments undergo collective vibrations because of in¬ 
ternal excitation. The collective excitation modes of the 
fragments may be investigated using both the time and 
the frequency domains. To perform Fourier analysis, we 
make use of the spectral power function 

P c (w) = [Re C(w)] 2 + [Im CM] 2 , (10) 

where C( w ) is the Fourier transform of the moment of in¬ 
terest. Unfortunately, owing to the limitations in the nu¬ 
merical grid size (the z direction spanned 160 grid points, 
ranging from -79.5 to 79.5 fm), a signal corresponding to 
the evolution of the multipole parameters of the individ¬ 
ual fragments could only be measured for approximately 
1000 fm/c before the grid boundaries were approached. 
For a signal of this length, the resolution of the calcu¬ 
lated power spectrum is of the order Hlo = ~ 0.6 

MeV [69]. 

To extend the time evolution domain and consequently 
improve energy resolution, we adopted the following ap- 
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proach. Rather than performing the calculations in an 
impractically large numerical grid, a Galilean transfor¬ 
mation is applied to each fission fragment. The boost 
momentum is chosen to cancel the corresponding linear 
momenta of each fragment. After boosting, the two frag¬ 
ments remain approximately still in the box, and their 
excitation modes can be studied for a much longer pe¬ 
riod. The evolution of the configuration with an ini¬ 
tial /? 2 o = 1-25 is presented as an example of this new 
method. 

Inside the masked regions of space corresponding to 
the fragments, the linear momentum may be calculated 
by integrating the current density: 

Pfrag = J dr j{r) ■ (11) 

This momentum therefore has units of velocity [30] . The 
linear momentum of the fragments may then be instan¬ 
taneously removed by applying a Galilean boost to the 
single-particle wave functions, 


<p(r) = exp ^ i P r j <p(r) , (12) 

where Af rag is the integrated particle density correspond¬ 
ing to the fragment and tp(r) are the single-particle wave 
functions. The Galilean transformation should be ap¬ 
plied in the masked region of space with the correspond¬ 
ing momentum for each fragment. The effect of the trans¬ 
formation is to effectively boost the particles in the oppo¬ 
site direction with the exact momentum they are propa¬ 
gating with through the grid. 

Figure 16 shows the decomposed energy functional for 
the fissioning case with initial /? 2 o = 1.25. The Galilean 
transform was applied when the separation of the frag¬ 
ments reached 100 fm (corresponding to f ~ 1900 fm/c), 
and the calculation terminated at separation 105 fm. 
Upon application of the transformation, panel (j), cor¬ 
responding to the total energy, displays a decrease of ap¬ 
proximately 140 MeV. This mirrors the drop in collective 
kinetic energy by the same amount, which corresponds 
to an instantaneous removal of excitation energy owing 
to translational motion. This demonstrates that the to¬ 
tal excitation energy is dominated by contributions from 
translational motion, rather than internal collective ex¬ 
citations. As for the other energy contributions, the nu¬ 
clear and Coulomb energies remain unaltered before and 
after the boost is applied around t « 1900 fm/c. This 
is to be expected, because these terms are all Galilean 
invariant. 

The collective kinetic energy drops instantaneously to 
ps 1.1 MeV following the transformation, as seen in inset 
panel (i). Because the translational energy has been re¬ 
moved at this point, the remaining collective energy can 
be interpreted as the sum of the internal excitation en¬ 
ergy shared between the two fragments. Reference [18] 
discusses an alternative method to deduce the internal 
collective excitation energy of the fragments, but the 
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FIG. 16. (Color online) Decomposed contributions to the 
energy of the system for the case with initial deformation 
/?20 = 1.25. A Galilean transformation has been applied to 
remove the linear momentum of the individual fragments once 
the separation between the fragment center of mass reaches 
100 fm. The calculation is terminated once the separation 
exceeds 105 fm. 


method applied assumes a priori knowledge of the fis¬ 
sion products of the system. 

The internal collective excitation energy is very small 
compared to the total excitation energy released in the 
fission process, which is ~ 180 MeV (see Fig. 14, /3 2 o = 
1.25). This justifies the previous assumption that the 
final collective excitation energy deduced in TDHF is 
dominantly translational kinetic energy, so it may there¬ 
fore be compared to the experimentally measured kinetic 
energies (Fig. 15). Further, the energy functional con¬ 
tributions (Fig. 16) may be compared to these in Fig. 
10 , where the calculation was terminated at the point 
where the transformation is applied in this case. Figure 
16 demonstrates that the nuclear potential part of the en¬ 
ergy functional is unaffected by the transformation. The 
calculation was performed in a grid of identical dimen¬ 
sions to those presented in Fig. 10 (48 x 48 x 160 points), 
and the time elapsed has effectively doubled from those 
previous calculations. As the measurement time of the 
post-fission fragments has been elongated, the resolution 
of the resulting power spectra will be enhanced accord¬ 
ingly. 

We note that the Coulomb interaction is long-ranged, 
so that even at a separation of 100 fm there is an in¬ 
teraction between the fragments. Translational motion 
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resumes after the Galilean transformation is applied, 
and the translational kinetic energy slowly increases. 
This can be seen by the gradual increase of the collec¬ 
tive kinetic energy in Fig. 16 following the transforma¬ 
tion. Therefore, the center-of-mass separation eventually 
reaches 105 fm and the calculation is terminated. We 
note that, in principle, one could reapply the Galilean 
transformation at every iteration to extend the total 
time and improve energy resolution. This artificially pro¬ 
longed Coulomb interaction between the fragments does 
not alter significantly the dynamics of the fragments. 
Any potential adverse effects have not been studied in 
detail, and alternate methods can also be applied to mit¬ 
igate this increased Coulomb interaction [70, 71]. 



FIG. 17. (Color online) Evolution of (a) quadrupole, (b) 
octupole and (c) hexadecapole deformation parameters for 
the heavy fission fragment. The initial deformation was 
P 20 = 1.25. The measurement time is significantly extended 
by applying the Galilean transformation to remove the linear 
momentum of the fragments. See text for more details. 

The evolution of the multipole deformation parame¬ 
ters for the two fragments are shown in Figs. 17 and 18, 
corresponding to the heavy and light fragment, respec¬ 
tively. The corresponding power spectra are shown in 
Figs. 19 and 20. Let us stress that the resolution of the 
spectra is ~ 0.4 MeV, a significant improvement with re¬ 
spect to what would be obtained without extending the 
measurement time. 

The evolution of the multipole fragments is qualita¬ 
tively similar for both fragments. The quadrupole os¬ 
cillations are centered around values of (3 2 o ~ 0.18 and 
0.35 for the heavy and light fragments, respectively. The 
corresponding hexadecapole moments oscillate in phase 
with the quadrupole deformation around central, non¬ 
zero values. The corresponding octupole deformations, 
in contrast, are modulated around a zero value, and it 



FIG. 18. (Color online) Evolution of (a) quadrupole, (b) oc¬ 
tupole and (c) hexadecapole deformation parameters for the 
light fission fragment. The initial deformation is /?20 = 1.25. 



FIG. 19. (Color online) Power spectra corresponding to Fig. 
17 (heavy fission fragment). The resolution is significantly 
improved due to the longer measurement time available with 
the use of Galilean transformations to remove the linear mo¬ 
mentum of the fragments. 


is difficult to ascribe a final octupole deformation for 
these (excited) fragments. While the specific deforma¬ 
tions might not be particularly relevant, the excitation 
patterns of the multipoles provide information on the 
collective vibrations of both fragments. We note that, 
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FIG. 20. (Color online) Power spectra corresponding to Fig. 
18 (light fission fragment). 


whereas a pattern of well defined, relatively rapid oscilla¬ 
tions are found for the light fragment, the heavy fragment 
multipoles have a more erratic time evolution. 

These features are reflected in the corresponding power 
spectra. We insist here that these can be obtained only 
with the required resolution after a boost has been per¬ 
formed on both fragments. Within the resulting spectra 
presented for the heavy fragment in Fig. 19, there is a 
well-defined peak for each multipole parameter between 
ss 1 and 3 MeV. For the light fragment, in contrast, Fig. 
20 shows a very well-defined peak for each multipole pa¬ 
rameter between 4 and 6 MeV. The latter reflects the very 
well-defined oscillations in the moments observed earlier. 
We speculate that these well-defined oscillations corre¬ 
spond to collective excitations within the two fragments. 
It would be of interest to compare the excitation modes 
obtained for the post-fission fragments to standard cal¬ 
culations of resonances in the corresponding nuclei. This 
would give access to a microscopic understanding of the 
fragment excitation properties, including their tempera¬ 
tures and phonon structure. 


V. CONCLUSIONS 

We have presented an analysis of the fission pro¬ 
cess using TDHF techniques as implemented in SKy3d. 
Starting from the calculated one-fragment quadrupole- 
constrained PES, the dynamics of fission were investi¬ 
gated. Deformation-induced fission was explored by re¬ 
leasing the quadrupole constraint and time-evolving a se¬ 
lection of states situated below, around and beyond the 
second static fission barrier. Three behaviours were ob¬ 


served. The states with a quadrupole deformation be¬ 
low the peak of the fission barrier undergo vibrations 
corresponding to a collective giant resonance. For these 
states, DIF is forbidden in TDHF, because a collective 
tunnelling through the barrier must occur to reach a fis¬ 
sioned configuration. 

A different behaviour is observed for the evolution of 
states which are situated beyond the peak of the sec¬ 
ond static fission barrier, but before the critical point 
where the static one and two-fragment pathways inter¬ 
cept. Upon time evolution up to 9000 fm/c, these states 
also fail to fission, but the dynamics are not typical of 
collective giant resonant modes. The repulsive Coulomb 
force attempts to drive the configuration towards a fission 
point, but owing to the competition with the attractive 
terms in the energy functional, scission does not occur. 
DIF is inhibited for these initial configurations, and it 
can only be speculated if these states would eventually 
fission with a longer time evolution. 

For states with a static deformation exceeding the in¬ 
tersection of the one and two-fragment fission pathways, 
DIF was observed upon time evolution. We interpret 
that, because a static two-fragment configuration exists 
with greater binding energy than the one-fragment con¬ 
figuration, it becomes energetically possible for the one- 
fragment configurations to evolve to fission with only a 
modest rearrangement of the densities. The evolution of 
the pre-fission fragment displays a rearrangement of the 
densities up until around the point of scission. At this 
point, the Coulomb repulsion between the pre-formed 
fragments overpowers the nuclear potential, and trans¬ 
lational motion sets in as the fission products rapidly 
accelerate away from one another. The timescale for 
DIF varies depending upon the deformation of the initial 
state. The least deformed configuration demonstrates a 
density rearrangement phase lasting approximately 1500 
fm/c, whereas the most deformed configuration is ini¬ 
tially close to the point of scission and the neck ruptures 
within 100 — 200 fm/c. 

A selection of fission products was observed for the 
various initial configurations considered. When com¬ 
pared to experimental measurements of neutron-induced 
fission processes, the agreement of the calculated frag¬ 
ment masses demonstrates promising results, although 
several effects, particularly dynamical pairing, are miss¬ 
ing in our approach. The energy released is shown to 
be dominantly translational kinetic energy, and agree¬ 
ment between theory and experiment was found to be 
reasonable when comparing the calculated and measured 
kinetic energies of the post-fissioned systems. We have 
pioneered a method to remove translational kinetic en¬ 
ergy of the post-scission fragments to provide details of 
their internal excitations. 

This initial investigation into fission induced by defor¬ 
mation effects using time-dependent techniques has pro¬ 
vided insightful results on the interplay between struc¬ 
ture and dynamics in the fission process. We have also 
devised a set of useful computational analysis tools for 
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the post-scission fragments. Deformation-induced fis¬ 
sion, however, provides a limited amount of two-fragment 
configurations. Further extensions of the present research 
involving particle-number projection could provide access 
to more relevant mass distributions. Moreover, a refined 
step-by-step linear momentum removal could help extend 
the time length and energy resolution of the two-fragment 
excitation data. 

The present analysis is relevant for a variety of fission 
processes. Spontaneous fission would presumably tunnel 
the system across the PES barrier into random states 
within the “allowed” region. These would subsequently 
decay into the two-fission pathway. One can in principle 


explore a wider landscape of fission fragments by exciting 
nuclei along their fission paths [50]. This theoretical ap¬ 
proach can be linked more naturally to induced fission, 
where the energy deposited by external probes induces 
the fission process. We plan to explore this approach 
using TDHF techniques in the near future. 
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